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The computer procian calculates tne supersonic f 1 c * ever 
a cere in air or water. The rain objective is to calculate 
tie cere semi-vertex ancle given prescribed initial condi- 
tions. Tie program is written in structured FCRTRAN and 
implements Busemann’s graphical integration technigue. 
Supersonic flow over a cone in water is useful as a good 
first approximation tc tne motion of tie metal jet frem an 
explosive shaped-char ce fired underwater. 

A topical result for supersonic flow over a cone in 
water is as follows: given an upstream temperature, 223.16 

Kelvin; upstream pressure, 1 tar; shock angle, 20.0 degrees; 
and pressure behind the shock front, 5 kiiotars, the cone 
semi-vertex angle is calculated to be 7.23 degrees. 

Generally, pressures involved in water flow are much 
larger than for air flow, and the cone semi-vertex angles 
for water flow are smaller than for air flew. 
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I. IN 15CEOC1IOH 



EESCBIfllON CP TEI PfiOBIEh 

lie solution of the hydrodynamic equations describing 
supersonic flow ever a cone in air has been well Known since 
the 1930's. Until recently, the problem of describing the 
flow ever a cone in water has been limited to solutions of 
the subsonic case. primarily, calculations were United to 
subsonic flew because researchers believed that surer sonic 
flow in water was not feasible for normal vessels (such as a 
ship). Per ordinary vessels in water, it is certainly true 
that supersonic flow past that vessel is highly improbable. 
However, the motion of the metal jet from an explosive 
shap e d-c bar ce fired underwater is supersonic. 

I his thesis presents a computer program which calculates 
the hydrodynamic flew past a cone in either water ci air 
under supersonic conditions. I he program utilizes the 
methods developed by previous researchers f or calculating 
the supersonic flow in air and which have been suitably 
modified to describe the conditions in the water. item 
modifications include utilizing the modified Tait eguation, 
which is the "thermal" or "thermodynamic" equation c: state 
for water, to describe the physical state of the water 
rather than the perfect gas law used for air. 

In actuality, the cone liner in the jet from an explo- 
sive sha p e a-cha rg e is blunt-nosed rather than conical. 
However, solution of the conical case is a preliminary 
requirement to solution of the actual blunt-nosed problem. 
Ihe scluticn to the conical flow case will serve as an 
excellent test program for the solution to the bluit-r.csed 
problem. Inis thesis presents a solution to the conical 
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c on tiruiig research into the problems or uuiiizing exp 
sha p e c-char ges ir ar underwater environment. 
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E. ailECDCICGY 

lie computer prccram presented in tnis work was origi- 
nally developed in the BASIC computer language using a 
Hewle tt-Fackar d HP-61 computer. That program is the tasis 
for this thesis. It was desired to translate the pregram 
into a higher-order computer language for executicn cn a 
large, mainframe computer system. Ihis translation was 
desired in order to make the program more readily accessible 
to a wider body of researchers and in order to speed the 
execution time of the pregram. In this thesis, the 
relieving goals have been accomplished: 

• Successfully translate the program frem BASIC into a 
higter-erder language. This goal was met by utilizing 
FCFTFAN as the high-level language of choice. T *hile 
FC5TFAN has many drawbacks as a high-order language, it 
is still widely used in the scientific community. 
FCFTFAN was used, therefore, so that the program will be 
useful to as wide an area of researchers as pcssihie. 

• Follow modern programming practices in the design and 
implementation of the program. As before, the choice cf 
FCFI.FAN as the high-level language maxes this goal some- 
what acre difficult. However, many computer scientists 
have demonstrated that structured programming practices 
can he achieved using FCBIBAN. To the largest extent 
possible structured programming practices have been 
utilized. 

• Present a "user-friendly", well-documented program. In 
this regard, literal use of comments occur in the 
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A. EACKGECOHD 

lie fiist step necessary tc describe supersonic flow- 
over a core is to calculate tie tier mod yna mic properties 
across a shcck wave. The basic research into the chance of 
tier mcdyranic properties across a shock wave in water was 
extensively conducted and reported upon in luc^rwater 
Explosion Re sea rch £Ee£. 1] during and gust after acrid rtar 
II. lie test s ummarization of tnese works can re fcurd in 
Cole £ Eef - 2], 

lie primary source used as reference for the calculation 
ci tie t ydr cdy na lie properties of sea water at the frcr.t of 
a shcck wave is the werk of Richardson, Arons, and Halverson 
£Ref. 3]. They utilized graphical techniques, which were 
rather crude and te cic us, tc calculate the thermcayranic 
data needed tc describe the conditions of the sea water. 
Euhs £Ref. 4] used tie work of Richardson, et_ai-_, £Eef. 3], 
to develop a computer program for the HP41CV hand-held 
calculator which efficiently calculated the same thermody- 
namic properties. funs’ £Eef. 4] programs provided the 
lasis for the FORTRAN subroutines, used in this work, which 
calculate the same tier mod y na trie properties such as pres- 
sure, temperature, arc density. 

The fundamental equations used to calculate the tnermo- 
dynamic charges which cccur across an obligue shock frert in 
air have been well krewn fer decades. These equations are 
excep t icra lly well described in Kinney and Graham £Eef. 5] 
and form the basis for the initial calculations fer the 
supersonic flow past a cone in air. 
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ramie state ci air is the perfect gas lav. Ir vater, ere 
modified Tait equation is the ecuatioc most often used to 
cescriie the thermcdyramic state of t hs water. The modified 
lait equation can be used tc describe either pure water or 
sea water. The fori of the modified lait equation used in 
this thesis was taxer from Holt [fief. 6 ] , who has continued 
to perform research ir underwater explosion phenomena. An 
excellent discussion of the modified Tait equation ar.d how 
it can be utilized is contained in fiowiinson [fief. 7*. 

After tie thermodynamic properties of the water (or air 
as the case may be) on the downstream side of the shock 
front ha-56 teen calculated, an iteration method, utilizing 
an automated graphical construction first developed by 
Eusemann, [Eef. 8], is used tc progress from the shock front 
to the cone surface. The equations needed for use by this 
iteraticr method are fully described, for air, by Shapiro 
[fief- 6] # who describes their development and use. The 
methods which apply tc the flew past a cone in air can, with 
the .necessary changes made for the differences in the tner- 
mody namics of the two fluids, be used to calculate the 
conical supersonic flow in water- These calculations form 
the basis for the main part of the FCBTfiAN program which 
follows. It should be noted here that Shapiro also pcirts 
cut tie pioneering work of Taylor and b accord in the 1S3C's 
and 1S40's on the methods of solution of the flew over a 
cone in air problem £Fef. 8j. 



£. GZNEEAI EQUATIONS 

Ir the development of the equations which follow, it is 
assumed that these equations can be validly used tc calcu- 
late the f low over the cone in either air or water. The 
equations were primarily developed by Shapiro [fief. 8[ for 
flow in air. 
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as a result of the shock front, it is assured that the fluid 
prooeiti.es will remain constant on imaginary cones having a 
ccmmcr vertex. By this assumption, the flow past the cone 
can he calculated. The flew geometry is illustrated by 
figure 2.1 ana figure 2.2, which show a typical streamline 
and its image m the hedo graph plane. 

As discussed thcicughly by Shapiro [ Eef . 8], there is a 
discc i t in uc us change in noth direction and velocity acrcss 
the cbiicue shock frert. Feints 1 and 2 (see figure 2. 1 and 
figure 2.2) lie, therefore, cn a common shock pclar which 
originates at point 1. Between point 2 and the point ’s’, 
which is cn the cone surface only at infinite distance, 
there is a regicn of conical flow where the stream proper- 
ties vary ccntinuously. The velocity vector to the point 
’s' in tie holograph plane defines what is called the ccne 
semi-vertex angle with the centerline axis. In the methods 
which fellow, the ccne semi-vertex angle is the variable 
which is ultimately determined. Further, as pointed cut by 
Shapiic [fief. 8], since all streamlines in the flew experi- 
ence the same entropy jump acrcss the shock front, the flew 
between tie shock ficnt and the cone surface is beth isen- 
tropic and irrotationai. 

The second-crdei differential eguation which actually 
describes the flow cf the fluid past the - cone in air is 
fully developed n y Shapiro [fief- 8]- Shapiro notes that 
there are two methods commonly used to solve this eguaticn. 
Cne method, developed first by Taylor and Maccoll [fief. 8], 
performs a rumerical integraticn of the eguation. The second 
methed, which integrates the eguation using a graphical 
construction method, was first developed by Busemarn. The 
program developed in this thesis utilizes Busemann's methed, 
modified fer performance on a modern- high-speed computer, to 
perform the integraticn of the full second-order differen- 
tial eguaticn. 
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Figure 2-2 Hodograph Image of Typical Streamline. 
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derivation cf the equations can be found in Shapiro 
[Ref. S], and need net be repeated nere. The n omenciature 
used in the equations developed in this cnapter is detailed 
in 3a tie I. In keeping >i i t h modern thought, the h-K-3 
(meter-kilcgram-secord) unit system has been used throughout 
this thesis except for occasional lapses during the water 
calculations when pxessures are referred to in units of 
kilci a is . 

ricm the geometry of figure 2.3 , it can be seen that; 

V = V ccs(o) - 9) and V = -V sin {* - 6) (2.1) 

r 03 

In the development ci the actual second-order differential 
equation, Shapiro [Refo 8] shews that, due to the condition 
cf iriotaticnality , the following relation must be true: 



dV 



r 
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lif f er entia ting E^n. 2.1 with respect to 



one c 1 1 a i n s ; 



-r^~ = — V ( 1 - — ) sin (oo - 0) + — cos(oo - 9) 
a 03 aOJ dco 



( 2 . 3 ) 
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figure 2.2 Cone and Flow Geometry. 
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dV 



= V = -V sin(^ - :) 



( 2 . n ) 



Substitution of this result into Egn 2.3 yields: 



dV 

dno 



dQ _ 

doo Vtan(co - 0) 



( 2 . 6 ) 



Shapiro [fief. 8] shots that the equation governing the flew 
in the ccrical region is: 



k - 1 dV 03 9 2 2 

( - --) (2V + V cotoo + -r-^)(V“ - V - V ) 

v 2 r oj dco max 03 r 

dV AV 

= (V -r-^ + V dV co 
r doo co — )V 

dco co 



(2.7) 
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Eliminating — from Zgn 2.4 and substituting the excres- 
duj dV r dV M 

sions for V ,V , -j— - and — 4 

r (o' dco doT 

2.5, into Egn 2.7 gives: 



given by r,gns 2. 1 , 2.4 and 
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sin6 

Vsin(oo - 0) sinoo 
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2V 2 sin 2 (M - 0) 

(k - D(v 2 - v 2 ) 

max 



(2.5) 



Eesignating 1 E ’ as tie radius of curvature of the hodccrap] 
streanlire, one obtains: 



R = 



dV 



sin(co - 0)dto 



( 2 . 5 ) 
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ccn-cmii 3 iqr 2.5 »itn Zqr 2.3, ore arrives at : 
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2V sin“(r - 9) 
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But, the energy equation asserts that: 




( 2 . 11 ) 
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where c is the local speed of sound in the fluid. 

Equation 2.12 is the basis for the calculation of the 
supersonic flow ever the solid cone in either air or water. 
The graphical integration method invented by Busemann is 
adequately explained in Shapiro [Ref. 3] and need net be 
repeated here. Essentially, the computer program given in 
this work automates the Eusemann graphical integration 
methcc for calculating the cone semi-vertex angle. 

C. ECOATICKS SPECIEIC TC AIR 

The equation of stare for air is specified by the 
perfect gas law (under the assumption, that is, that the air 
hehaves as a perfect gas). This law is quite elegant and 
simple ana allows eas^ manipulation to obtain various quan- 
tities. The form cf the perfect gas law mest often used in 
the calculations of this thesis is: 



p = pET 



(2. 13) 
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(2. 14) 



In equation 2.14, M is the molecular weight ox the air. 

Ir air, the charge in the thermodynamic properties of 
the cas as it cresses the shock front are easily calculated, 
as shewn in Kinney aid Graham £fief. 5]. Since most fluid 
dynamics textbooks illustrate the develo r ment or the equa- 
tions which follow, it is net necessary to derive them here. 
3s mentioned previously, Kinney and Graham [Eef. 5] provide 
exceptionally lucid explanations and derivations. The prin- 
ciple equations used to calculate the thermodynamic chances 
which occur across the shock front in air are as follows; 



tan(S - 6 ) 
f tanS 



2 + (k - 1) M“ sin“3 
(k + 1) sin“B 



<2. 15) 



Eauaticn 2.15 is used to iteratively determine the deflec- 
tion ancle 9 . All ether quantities in this equation axe 
known (i.e. 3 is the snoex angle and M ^ is the freestream 
Mach number, both cf which are input parameters to the 
urogram for the air c alcula tic ns) . 

2 + (k - 1) M 2 sin 2 3 

[M 9 sin(3 - 9)] = y 5 (2. 1 c) 

2kM^ sin 3 - (k - 1) 

Having determined 9 from eguaticn 2.15, equation 2.16 is 
used tc determine the Mach number on the downstream side of 
the shock front. 
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figure 2.4 Gecietry cf the Obligue Shock Front. 

2 2 



I. ZCGATIC.NS 3PZCIE1C IC 5 ATZP 

Ir contrast to air , which has a relatively elegant a nd 
simple state e-uarici, the equation or state for water is 
rather acre complicated. The nest commonly used state equa- 
tion for water is known as the modified Tait equation, which 
may re written as: 

P = B (S ) [ (■=•) n _ i] 21 ) 

P 



where £ (Z) is a slowly varying function of entropy aicne, n 
is ap qr cxim ately a ccrstant equal to 7.15, and c is the 
value cf the density for zero pressure. The above equation 
is from Eclt [Eef. 6], but, in different forms, it is also 
described quite extensively in Cole [Ref. 2], Richardson, 

et^ a 1_. , [Eef- 3] ard Eowlinscn [Eef. 7]. As mentioned in 

the Introduction to this work, the modified Tait equation is 
the "thermal’' eguaticn of state for water. Cole [Eef. 2] 
shows that the modified Tait equation is of the form 



P 



B (S) 



[( 



v (T , 0) v 
V (T , p) 



n 



1 ] 



( 2 . 22 ) 



or, in simpler terms. 



p = p ( v , T ) 



( 2 . 23 ) 



Since the full modified Tait equation relates the three 
ther medy na a ic quantities of p, v, and T, it is called the 
"thermodynamic" cr "thermal" equation of state. 



Similarly, ir ccrcrast tc air, tne thermodynamic cranges 
which cccur as a result of the shock process in water cannot 
he easily delineated iy simple equations as in the air case. 
However, a s impiif ica tion can he made in the water case 
because, urlike in the air case, the pressure jump across 
the shock in water is so very large. In air, pressure 
changes across the shcck front on the order of 1 tc 2 bars 
are considered large (at least for chemical explosions) . In 

contrast, as pointed cut by Richardson, et. al^, [Eef. 3], 

the pressure jump across a shcck in water is on the order of 
kilctars tc tens of kiiohars. Therefore, the calculations 
can he simplified by specifying the pressure on the down- 
stream side of the shoe k front. This is valid since the 

upstream pressure is sc small in comparison to the upstream 

2 

dynamic pressure 'V 1 /2 and in comparison to beth tne 
downstream pressure and dynamic pressure. The specification 
cf the downstream pressure is accomplished, in the program 
of this thesis, by allowing the user to input a "pressure 
multiplication factor {MfACT)", which converts the pressure 
upstream tc a pressure downstream at point 2 which is grven 
by: 



p = p I ME ACT X 1010.0 (2.24) 

t 2 n 

The factor 1000.0 ir eguaticn 2.24 converts tne richt-hand 
side cf the eguaticn from a pressure in bars to a pressure 
in kilotars. As an example, if p ^ is 1 bar and MFACT is 5.0, 
the pressure at point 2 downstream will be 5 kilobars. 

It can be shown that the simplifying assumption made 
above is entirely valid by considering the momentum eguaticn 
for steady frictionless flow along a streamline. The 
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cre-dunersicnai momentum equation stares: 
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d + c V - / 7 = d + 0 V - / 7 
^ 1 l i'“ • 7 2 7 ' “ 



Now r:ct€ that tjpical values for and V for tie fret lets 
under consideration are as fellows; 

P 1 = 1 000 kg/m and V = 1500 m/s 

Therefore, /2 = 1.12 X 1 £j Pascals or about 11.1 kilo- 

tars. Tjpical values for p , the upstream pressure, are on 
the order cf 1 to perhaps a few tens of tars. Thus, it can 
re safely assumed that p can he ignored compared to tne 
upstream dynamic pressure. 

Having specified this pressure, the program cf tr.is 
thesis utilizes the subroutines developed by Fuhs [Ref. 4] 
to calculate the density and velocity components at a point 
just downstream cf the shock front. Since the geometry of 
the velocity flew across the shock front is the same in 
either air cr water (see figure 2.4) , the velocity compo- 
nents determined in the first step can be used to "jump back 
across" the shock front to the upstream side where the free- 

stream Mach number M can be determined. 

1 

After determining the conditions on the downstream side 
cf the shock front, the program iterates along the stream- 
lines to determine the cone semi-vertex angle. This itera- 
tion is performed using the Busemann graphical construction 
as for tie air case. The major difference between the water 
and air cases, after the initial shocked conditions have 
teen determined, is that the local speed cf sound on a 
streamline must be determined iteratively, for the water 
case, bj using Fuhs’ [Kef. 4] subroutines to determine the 
density and pressure at each point. Knowing the pressure 
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and density, the speed of sound m to a water can re calcu- 
lated ire a : 



= (^)v4) n 



(2 . 26 ) 



which is derived free the modified Tait equation using the 
definition cf the speed of sound as: 

2 _ 9p_i (2-27) 

c 3p ' S 



Ihese equations, combined with the equations developed in 
secticn E cf this chapter, constitute ail the equations 
needed tc completely solve the supersonic flow over a ccne. 
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r Sacius m spencal cccraiaa tes 

f Radius cf cavature of the 

hcdograph streamline 

1 lemperature 

V Velocity 

V Maximum velocity for adiabatic flow 

max 

3 Angle the streamline maxes with the shock 

plane (same as the angle c ) in air caiculaticns 

6 Ccne semi-vertex ancle 

9 flew direction 

p Mass density 

o Shock angle 

a) Ancle in spherical coordinates 
A Universal gas constant 

( ) Signifies a condition upstream of 

1 conical shock or a sheex front 

( ) Signifies a condition downstream cf 

2 conical shock cr a shock front 

( ) Signifies conditions at cone surface 
( ) Signifies a component in r-airection 

( ) Signifies a component in oj-direction 

oo c 



TABU I 
Nom encla t ure 



c Speed of sound 

k Ratio of specific heats (c /c ) 

p V 

M Mach number 

M Molecular weight of air 

a 

p Pressure 
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I II. DESCRIPTION CF THE COMPUTES PROGRAM 
A. PRCGEAb LOGIC 

As discussed in Chapter 1, the computer program devel- 
oped ici this thesis »as written in the FORTRAN programming 
language. Structured programming practices have teen 
followed throughout in that the program is divided into 
blocks (or modules) of code, each of which is designed to 
perform a single calculation sequence. Many or the modules 
have leer placed inline rather than being writter as sepa- 
rate subroutines or functions. This was dene prinarilg for 
ease ci use and understanding. In addition, however, most 
modules were placed inline because the parameter list 
exchanged between the module and the main pro gram would have 
teen excessively long otherwise. Regardless of whether the 
code is included inline or as a separate module, the program 
was writter in a manner which ensured that the "side- 
effects' 1 problem discussed by MacLennan [Ref. 9] did ret 
cccur. In addition, the use of FORTRAN'S COMMON construct 
has been studiously avoided to prevent the aliasing problem 
discussed by MacLennar [Ref. S ]. A module hierarchy chart, 
which shews the mayor modules of the program and their 
interconnections, is given as figure 3.1. 

The flew logic of the main computer program is illus- 
trated ly figures A. 1 through A. 7 of Appendix A. These flow- 
charts shew that the initial part of the program (up to line 
60) is used to initialize certain key variables and to 
gather input values from the program user. Note that, as 
stated in the objectives for this program, all user input is 
verified tc ensure that the values entered fall within 
prescribed limits. In erder to ensure that input values are 
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Figure 3.1 Module Hierarchy. 




witnir limits, the surroutine CHKINP is irvoKed and the 
input values are passed as parameters. Ihe flew lcgic of 
subroutine CHKINP, whicn is se 1 r-ex plana tor y , is ircluded as 
figures A. 8 through 1.10 of Appendix A. 

After receiving and verifying the user’s input, the 
program begins the calculations required to determine the 
cone semi-vertex ancle. from line 60 to line 40C, the 
program calculates the initial thermodynamic properties of 
the water cr air, as the case may be, and calculates the 
changes which occur in the fluid properties as a result of 
the passage of the shock front. As mentioned previously, 
these calculations fer the air case are ratner straightfor- 
ward and, because of this simplicity, the computer code for 
the air calculation has been placed inline. For the water 
case, the subroutines WSHOCK and EPV are needed to calculate 
the required thermodynamic properties. Subroutine WSHCCK 
calculates the ther mcdynamic properties of the water at 
point 1 upstream of the shock front and at point 2, just 
across the shock frcr.t on the downstream or shocked side. 
Appendix A, figures A. 12 and A. 13, contains the flowchart 
which describes the logic of this subroutine. Subroutine EPV 
was copied, with the author's permission, from Fuhs £Eef. 4] 
and is fully described in that work. However, the ECEIHAN 
translaticn of EPV is included in the program listing 
contained in Appendix C. 

Having calculated the initial thermodynamic fluid condi- 
tions cn both sides of the shock front, the program begins 
the iteraticn process required to determine the cone semi- 
vertex angle. These calculations begin at line 400 of the 
computer program. This iteraticn process makes use of the 
Eusemann graphical integration technique discussed by 
Shapiro £Eef. 8J. Eirst, the radius of the 3usemarn curve 

for a pcint J is calculated using the velocity at that 
point, the streamline angle at that point, and the value cf 
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tne angle omega at that point (see figure 5.2). Next, tie 
center of the circle for the Eusenann curve at the point J 
is calculated using the velocity con toner, ts at point J 
(which are designated u and V) , the radius calculated ir the 
previous step, and the value cf the angle omega. This circle 
center provides the pcint from which the program "draws" the 
arc used to calculate the next point on the Eusemann curve, 
rinally, tie velocity and the Eusemann velocity components, 
at the next point, J*1, are calculated by "swinging an arc" 
from the circle center calculated in the previous step. In 
addition, the streamline angle at the next point is 
calc ula t ea . 

The program next tests to see if the cone surface has 
teen reached (line 475) . This test is conducted by deter- 
mining if the absolute value cf the difference between the 
streamline angle and the angle omega is less than a speci- 
fied test value (which is set to 1 X 10 ~^) . If this differ- 
ence is less than the test value, the cone surface has been 
reached. In this case, the program then calculates the ther- 
modynamic fluid conditions at the cone surface and the ccne 
semi-vertex angle arc displays the final results cf the 
program. If the dif f erence between the streamline angle and 
the angle cmega is creater than the test value, the c one 
surface has not beer reached. In this case, the program 
calculates the th er mcay nam ic fluid properties at the next 

point, J+1, then loops bach to line 400 to begin the 
Eusemann graphical integration process for the new point. In 
the case of a water run, the subroutine WAIVE! is invoked 
just pricr to looping bach to line 400. This subroutine 
calculates the t hermcaynami c properties of the water at any 
given point. Appendix A, figures A. 14 and A. 15, gives a 
flowchart which demonstrates the logic flew witnin this 
sub routine. 
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Ficure 3„ 2 Graphical Construction cf Cone Flow 



The ficcrai includes a feature which allows the user to 
rake repeated executions of the program without the reguiie- 
ment cf "reloading" the program between each execution, This 
feature was included i y soliciting a response from the user 
at lire 15C0 as tc whether he/she wishes to sake another 
execution cf the program. An affirmative response causes the 
program tc loop tack to statement 1 at the beginning cf the 
Initialization section of the main program, A negative 
respcrse causes program ter mir a ti on . 

Ncte that certain separately compiled functions included 
with this program have not teen flowcharted because their 
logic is sc str aig h tf or war d and because the subprograms in 
guesticn usually consist of only one or two lines of execu- 
table code. The functions fitting into this category are: 
ETOfi, RICE, and 50SA1E. 

E. 0 SEE INS IE UC 1I0ES 

As discussed in the objectives listed in chapter 1, one 
goal cf this thesis was to ensure that the program presented 
was easy tc use and maintain. Ease of use has been facili- 
tated by including code which verifies that the input 
provided by the user is within specified ranges. Eor 
example, in the water calculations, the pressure at pcint 2 
can he nc .greater than 100-0 kilcbars. This is due tc the 
fact that funs’ [Ref. 4] subroutines, which are used to 
calculate the t he r m caynami c properties of the water, are 

based cn the work of fichardson, et . ah, [Ref. 3], which 

only gives results up to pressures of 100.0 kilctars. 
Therefore,, while the subroutines could compute the wiater 
conditions at pressures greater tnan 100.0 kilobars, it is 
uncertain whether the values so calculated would be entirely 
correct. I cr this reason, the range of the input values has 
been iestricted. Ease of use is furtner facilitated by 



pro v i c i n g clear, meaningful cut:at. Appendix 5 contains 
sample outputs generated by the r rogram for the suciar] and 
complete print cpticrs of the program fcr horn water ana 
air . 

Cce ci the irajoi criticisms leveled ry computer scien- 
tists against the FORTRAN language is fcr its lack cf a 
requirement for formally defining all variables used within 
a program. Similarly, FORTRAN has an implicit declaration 
policy thereby any variable whose first letter is I, J , K, 
I, E cr N is implicitly declared to be of type Integer. In 
this pregram, these two unfortunate char act er i s t ic s of 
FORTRAN are avoided by requiring explicit declaration or ail 
program variables. 

Ease of mairtenarce and ease of understanding cf the 
program have been achieved through the use of a literally 
commented program and through the use of relatively mean- 
ingful variable names. In this regard, the FORTRAN language 
is less than desirable since it limits variable names tc six 
characters. Appendix C contains a fully documented listing 
cf the computer program and all subroutines or functions 
used by the main program (ether than standard library .func- 
tions such as sir) . Tables VI, VII, VIII and IX, located in 
Appendix A, contain lists cf all variable names used in tne 
main program, their meaning, and their MKS (met er-kilcgrac- 
sedena) units, if appropriate. Table X contains a similar 
variable list for the subroutine KSHOCK and Table XI 
contains a variable list fcr the subroutine WATVFI. These 
two tables are also lecated in Appendix A. 

Finally, this pregram was developed using the FORTRAN IV 
language supplied by the IEM Corporation as part cf their 
IBM 3VC/3033AP computer system which is the main computer 
system available at the Naval Postgraduate School. As far as 
is krewr, no implementation specific features have been 
included in this program. Therefore, the program should 
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execute cl ary system which has a FORTRAN compiler. In crcer 
tc execute the program on the 153 370 computer systei avail- 
able at the Naval Postgraduate School, the rollcvirg steps 
rust le acccrplished it the order given; 

(1) lie program must te compiled using the conanc: 

FORTH 1 CCNEFiCW 

This command invckes the FORTRAN H Extended compiler, which 
is an optimizing, production run compiler supplied as part 
or the IBM computer system. The program could also he 
compiled using another FORTRAN compiler such as the FCRTZGI 
compiler. Note that the above assumes the program supplied 
by this thesis has teen entered into a file which has the 
filename CCKEFlOft and which is of filetype FORTRAN. Note 
further that this step need only be performed ere time 
provided errors do net occur. 

(2) Next, the libraries of standard subroutines and 
functions supplied by the computer center must be attached 
to the file. This is accomplished by issuing the following 
com m a r c : 

G 10 5 A I TXTLIZ FCRTMOD2 MGD2ZZH 

Note that this step need only be performed once provided 
abnormal terminations do not occur. 

(3) In the FORTRAN language, a formatted input/cutput 
(I/O) statement requires the definition of an I/O device. 
For example, at most installations, by default, I/O device 5 
is used for input and I/O device 6 is used for output. To 
utili2€ the program cf this thesis properly, the following 
commands must be issued to define the I/O devices tc be 
use d ; 
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FILIIrr 06 Z 12 : ( C C 'J r, F L G /i C j T ? U ^ 

FILEIZf 07 T Z R h (PER:!) 

As a result of the sieve commands, all cutout frci ire 
program will be written into a file on the user's disk named 
CONEfICii with a filetype cf OUTPUT. It is not strictly 
necessary tc issue tie FI1EEEF 06 command shown above, but, 
if net, all output will be written at the terminal and will 
not be available for printing. The second file defir.ition 

commard (IIIEDEF 07) is absolutely req uir ed or the program 

will ret operate. This file definition tells the computer 
that I/O unit 07 is the computer terminal. In the program 
of this thesis, all input from the user is reguested from 

I/O unit C7. Therefore, if unit 07 is not defined, the 

program cannot receive any input. 

(4) Finally, the program can be executed by issuing 
the command: 

ICAE COKEEIOW (START 

fchen procram execution is completed, a cryptic message cf 
the f cm R; T=C.01/C.01 1 6: 45: 00 will be displayed at the 
terminal. This is simply a message from the computer's oper- 
ating system indicating that the task just reguested has 

been completed. The user can now utilize the system’s oper- 

ating commands to review the output from the execution of 
the program . 

If the user desires to execute the program again, return 
to step (3) of the above procedure. If the output from the 
previous execution is no longer needed, the same FI1ZEZF 06 
command can be issued. Howrever, if the previous execution's 
output is needed, the f ilety pe of the FILEDEF 06 command 
should he changed (e.g. change OUTPUT in the above command 
to CCTPCI2) . 
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IV. PRCGBAtl CALCULATION RESULTS 



A. TECGEAh RESULTS ICE THE CALCULATION 01 CCNICAI flCS IN 

A IE 

1 1 € program accuxately calculates the flow over a ccr.e 
in air. Numerous program executions were made for the air 
case at various upstream Mach numbers and for various shcck 
angles. The results obtained from these program executions 
are summarized in Tables II and III. The numbers produced by 
tne program were compared for the variables shown in the 
tables to the tables given in Kinney and Graham [Eef. 5], 
against the graphs civen by 5hapiro £Bef. 8], and to the 
tables produced by Nepal [Eef. 10], The values given in 
these sources are included in Tables II and III, where 
appropriate, along with the calculated results. As can be 
seen from these two tables, the calculations made by the 
program cive ^uite accurate results (less than a 13 error in 
most cases). Certainly, the program is much more accurate 
than cne’s ability tc read the graphs presented in Shapiro 
[Eef. £ j . 

rased cn these ccnpariscns, it is believed that the part 
cf the program which calculates supersonic rlow over a ccne 
in air is accurate. The importance of this fact is that, by 
feeling confident that the procedure followed for the air 
case is valid and that this procedure has been correctly 
implemented in the programming language, it is safe to 
assume that the Eusemann calculation procedure utilized for 
the air case can be accurately applied to the water case 
provided the water conditions at each point are calculated 
cor r e c tly . 
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(3) The fcdlowmy variables were held constant at the values given for 
all r icy ram executions: T, = 298.16 Kelvin; P, = 1 01 300.0 Pascals; 
and ft, = 3.0. 1 1 



E. RFCGEAM RESULTS ECS THE CALCULATION OF CCNICAl FLC* IN 

SALES 

As mentioned in ths previous section, the calculation 
procedure was verified for the air case by comparison witn 
known results. This comparison showed the met nod used ir. the 
calculations was correct. Ihe next step taken was to deter- 
mine if the subroutines used to calculate the t n er me ay r.a mic 
properties of the water gave results comparable to these of 

Richardson, et. al, [Ref. 3], ana Fuhs [Ref. 4]. 

Calculations were performed which gave results which could 
be compared to Table II of Fuhs [Ref. 4]. Table IV presents 
a comparison between the numeric results calculated by the 
program and those given by Funs. As can be seen, these 
results are, within reasonable error, remarkably similar. 
This is net unexpected since the subroutines used by the 
main program to calculate these variables are essentially 
direct language translations of F'uhs 1 programs. Therefore, 
the results should be similar. 

Having verified that the procedure utilized in the 
calculation of the cone semi-vertex angle was ccrrect 
(through the air results) and having verified that the 
program correctly calculated the thermodynamic r rcperties cf 
the water at any point, it was believed that the program 
could be executed fer the water case, for various initial 
conditions, with certainty that the results so calculated 
would be accurate. Hcwever, a question arose as tc whether 
the calculations performed in the water case were indepen- 
dent cf the upstream thermodynamic properties of the water. 
In air, tie only quantities needed to calculate the ccne 
semi-vertex angle are the upstream Mach number ( M ^ ) and the 
shock angle ( a ) (e.g. see figure 17.7(a) of Shapiro 

[Ref. 8] ). Thus, for air, the calculation of the ccne 

semi-vertex angle is independent of the upstream temperature 
or pressure, and, therefore, universal curves can be drawn 
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Properties of Sea Water at a Shock Front 
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under all 



which descrioe the cone flew under all conditions a* is 
shown ty Shapiro [Bef. 8]. whether universal results nor 
the water case could he obtained was not sc clean; this 
point is new discussed. 

lueptrarn and Rcshho £Bef. 11] discuss the general ther- 
modynamics of fluids and introduce rhe concept cf a 
"theiflal" ecuaticn cf state. In general, the thermal equa- 
tion cf state fox any fluid is given by; 



f (P/P /l) = 0 



(4.1) 



for a perfect gas, the thermal equation of state is given 
ty: 

p v = El (4.2) 



where B = A /M . Another equation which relates tne thermody- 
namic variables e, v, and 1 is the "caloric” equation cf 
state, which, in a general form, is given ty; 

f { e , v , 1) = 0 (4.3) 



A 



perfect gas is defined by; 



e = c 1 

v 

cr h y ; 



(4.4) 



h = c 1 (4.5) 

where e is internal energy and h is enthalpy. Thus, the heat 
capacity (either c^ or c^ ) is a constant for a calcrically 
perfect cas. 



an: 



?. cihko [Bex- 11] dcicnsaate that 



tne 



i. d € j n a n r. 

thermal cic caloric €cuatioi£ cf state are tnermod 
related ty: 



amically 



3h ' 
3p 'T 



V 



T — I 

T 3T 'p 



(4 .b) 



Introducing equation 4.2 into equation 4.6 results in the 
right-hand side of equation 4.6 becoming e^ual to rero. 
Thus, for a thermally perfect gas, the enthalpy is indepen- 
dent cf the pressure, and, hence, enthalpy is a function 
only cf temperature. Therefore, a necessary conditicr. for 
equaticn 4.5 to be valid is to nave a thermally perfect gas. 

for water, Richardson, e t . a l . , £Bef. 3], demonstrate 

that the heat capacities are functions of the temperature cf 
the water. Thus, equations 4.4 and 4.5 are not valid for 
water, and water is ret a caloricaiiy perfect fluid. 

If a fluid is both thermally and caloricaiiy perfect, 
the conditions across a normal shock front depend only on k 
= c p /c v and on the freestream Mach number. Conversely, if 
the fluid is not thermally ard caloricaiiy perfect, addi- 
tional variables must be specified in order to define the 
shock conditions. Extending the argument, one can state mat 
the supersonic flow cf a thermally and caloricaiiy perfect 
fluid over a cone is a function only cf the heat capacity 
ratio, k ; the freestream Each number, ; and the shock 
front angle, q . As a result, a single graph is sufficient to 
represent all flews ever the cone. This is the condition for 
a perfect gas. 

In contrast, when the fluid is not thermally and calori- 
cally perfect, the solution for supersonic conical flew is 
dependent upon variables ether than k, M ^ , and c . 
Consequently , a universal graph- of the shock angle versus 
the freestream Mach rumber, with the cone angle as an entry 
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parameter, does not exj .so. Since water is net a thermally 
and calcricaliy perfect fluid, it is not possible tc plot 
universal graphs of tie variables of interest. 

To gain insight into the sensitivity of the supersonic 
conical flew in water to the upstream conditions, a study 
was conducted in which calculations were performed tc deter- 
mine whether the cone semi- vertex angle is relatively inde- 
pendent of the upstream temperature or whether various 
program executions are required to snow the variation of 
cone semi-vertex angle with water temperature. 

The ceterminatior discussed above was made by holding 
the upstream pressure, the shock angle, and the upstream 

Each number constant. Then the upstream temperature was 

varied tc determine the variation in the cone semi-vertex 
angle. In order to held Mp constant, the pressure multipli- 
cation factor, described in Chapter II. D, was varied until 
the Each number upstream for the new temperature matched 
(withir reasonable accuracy} the upstream Mach number for 
the original temperature. Note that "reasonable accuracy" 
for matching of the upstream Mach numbers meant matching the 
numbers tc the fourth decimal place. Table V provides the 
results from these c alculat ions . As can be seen from this 
table, the cone semi-vertex angles calculated for the 

various temperatures are all of approximately the same value 
(within 15»). The pressure downstream of the shock wave 
varies by 13% as can be seen in Table V. While there is small 
variation between the values for the cone semi-vertex 
angles, it must be remembered that, for prac tica l purposes, 
these variations are very small and can be safely ignored in 
examining the flow ever a cone in water. Thus, it was 
concluded that the calculation of the cone semi-vertex angle 
in water can be conducted independently of the water temper- 
ature (at least for the accuracy reguired for the calcula- 
tion of the movement cf the metal jet from a shaped-charg e) . 
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TABLE V 

Cone Semi-vertex Angle Variation with Temperature 
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litis ccr.clusion is true provides tiae upstream Mach rusher is 
held constant by varying the pressure rebind the shock ircnt 
tnrouch the mechanist cf the pressure multiplication factor. 
Kavinc determined that the calculations can he conducted 
independently of the temperature, plots of the variables of 
interest in water were cade in a manner similar tc the 
graphs for air shown in Shapiro [Ref. 8J. 

finally, having determined that the calculations could 
he cade independently of the upstream temperature and that 
these calculations wculd be accurate (within 1% fcr ccne 
angle and 11% for pressure) , it was decided to plot various 
paraceteis cf interest for the flow over a cone in water. 
Ihese graphs are presented as fellows: 

(1) figure 4.1 is a plot cf the freestream Mach Number 
versus the shock angle with the cone semi-vertex 
angle as an entry parameter. This graph is the 

water analog to Figure 17.7. (a) of Shapiro [Eef. £]. 

(1) figure 4.2 is a plot cf the freestream Mach runher 
versus the Mach number at the cone surface with the 
ccne semi-vertex angle as an entry parameter. This 
graph is the water analog to Figure 17.7. (c) of 

Shapiro [Ref . 8 jo 

(1) Finally, Figure 4.3 is a plot of the freestream Mach 
rumher versus the drag coefficient (c ) with the 
ccne semi-vertex angle as an entry parameter. This 
graph is the water analog to Figure 17.7. (f) of 

Shapiro [Ref. 8 j . 

It should be ncted here that the graphs developed and 
presented in this thesis were obtained through repeated 
executicns cf the C0NIF10W program, which generated approxi- 
mately 36GG data points per graph. These_ data points were 
pre- p r cc essed _ • t c remove entries with a Mach number of 
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greater tran 5.0 (a convenience :c: case cf presentation 

cniy) . 1h€ are- processed data points were passed tc an 

interpciaticn and smoothing program developed t y Franks 
£Bef. which used a guadratic Snepard's method tc sncctr 

the data feints generated. After smoothing, the data was 
post-processed to remove erroneous interpciaticn points 
which were generated because the original data was irreg- 
ular. Finally, the processed data was passed to a FCFTPA.N 
program developed ly the author of this thesis which 
utilized the Cortcuring feature of the DI5SP1A graphing 
routines of the ISSCC company. The results of this extensive 
processing are the giaphs presented in this thesis. Because 
cf the irregularity of the original data (which causes 
difficulties in the guadratic Shepard's interrelation 
method}, and because of the smoothing reguirea ty the 
contouring routines cf DISSP1A, the graphs presented are 
somewhat rcugh and should, therefore, only be used for "rack 
cf the envelope" calculations. Accurate results are 

provided ty the CONFFLOW program which should be utilized 
for mere precise work. 

Finally, one of the primary reasons for translating the 
original Easic programs into a higher-order language (in 
addition tc the desire to make the programs more accessible) 
was tc speed the execution cf the program. In this regard, 
the work cf this thesis has more tnan accomplished this 
result. It must be mentioned here that the program will 
execute extremely last, considering all the iterations 
needed, provided the computer system in use is net heavily 
in use at the time. For example, at the Naval Postgraduate 
School, when the computer system is relatively free, the 
program operates so cuickly that one execution reguiles less 
than 5 seconds (from the time cf the last user input tc the 
time cf the reguest by the pregram for an indication of 
whether another execution is desired). Conversely, when the 
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computer system is being utilized at near its capacity, the 
sane response tine car increase up to about 2 to 3 linutes. 
Ibis result, however, is sy at tomatic or any computer system 
being utilized at near capacity and is not a result cr a 
reflection of the program's structure or design. 
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figure 4.1 Freestrean Mach Number vs Shock Angle fcr Supersonic Flow in Hater. 
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7. C C NCI 0 £10 NS 



As discussed ir the Introduction, this thesis rad 
certair goals which it was desired to achieve. It is 
believed that these coals have all been successfully met. 
Ihe program of this thesis has teen successfully translated 
from Easic into FORTRAN as desired. In the translation of 
the original program, modern structured programming prac- 
tices have teen followed tc the greatest extent possible, 
finally, the program presented is quite "user-friendly” and 
is very well-documented. 

As was aem onstra t ed in the results chapter, this program 
calculates correct results for each of the air or water 
cases. It is mentioned here that, as pointed out frequently 
ty computer scientists, it is virtually impossible tc test a 
program for all possible cases. Therefore, no program is 
entirely error-free. It is believed that the program cl this 
thesis is as free of errors as is possible without exhaus- 
tive, time-consuming and extremely expensive testing. 

As a result of having met the three goals mentioned in 
the first paragraph, it is believed that the program of this 
thesis will provide an excellent working tcci for 
researchers in the field of shaped-charge jet penetration in 
water. Ihe fact that the program is well- doc u me n t e d will 
make any modifications or extensions to the program, should 
that be reguired or desired, easy to perform. Further, since 
the program is easy to use, only a cursory knowledge of 
computers is needed in order to utilize the program. These 
features are requirements for any computer program which is 
to be used by scientists as_a_tg_ol. Too often, computers and 
their programs require that the people who simply desire to 
use their capabilities must learn a great deal cf detail 



a tout tie computer order cc use toe cool. Is curs 
program, ai understanding of FCFIRAN is net reguired in 
order tc utilize the program cr its design. Unfortunately, 
due tc the operating system of the IB.i 370 computer system, 
the same cannot he said about the steps reguired in order to 
actually use the program. Seme familiarity with the computer 
system in use at the users location will he a necessity in 
order tc operate the program correctly. 

Finally, as mentioned in the Introduction tc this 
thesis, the program till provide an excellent test case and 
comparison model for a computer program which models the 
actual flew of the metal jet from an explosive shaped-charge 
fired through the water. Ihe actual situation is a hlunt- 
nosed, rather than a shar p-p cin ted conical, flew problem 
which is much more difficult tc solve. Therefore, a known 
solution and methodology of solution for the easier problem 
is a necessary first step to the solution of the larger 
problem. It is believed that the program of this thesis will 
serve as this necessary first step. Further, it is believed 
that, when the eguatiens solving the actual flow problem are 
developed in their final form, the program of this thesis, 
due tc its ease of modification, will serve as the program- 
ming model for the program which calculates the blunt-nosed 
flow problem. 
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This appendix certains the logic flowcharts of the main 
program and irs subroutines {i.e. those which have ret been 
described elsewhere or which are not part of standard 
computer center libraries) . The language is kept rather 
general sc the cverail logic cf the program can re demon- 
strated. If the reader desires to know how a particular 
logic seguer.ee is implemented, he need only refer directly 
to the program segment which the logical flowchart is 
describing. The flowcharts included in this appendix are 
listed below: 

(1) Ihe Main Program Flowchart consists of Figures A. 1 
through A. 7 . 

(2) Function CHKJhP consists of Figures A. 8 
through A. 1 0 . 

(3) Subroutine DIFANG consists of Figure A. 11. 

(4) Subroutine WSKCCK consists of Figures A. 12 and A. 13. 

(5) Subroutine WAIVE! consists of Figures A. 14 arc A- 15. 
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Figure A.1 Main Program Flowchart. 
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Figure A, 2 Main Program Flovchart (ccnt'd 
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Figure A„3 Main Program Flowchart (cent* 
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Figure A .. 4 Main Program Flowchart (ccnt'd 
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Figure A. 5 Main Program Flowchart (cent 
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Figure A. 6 Main Program Flowchart (cent 1 d. 
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Figure A . 7 Main Program Flowchart (cent 



1 



r 

[ 



I 



Ui co 

i 3 



«ll 




1 


*— 1 






rH 


3 


•H| 


1 


1 g 1 


I H 


1 G| 




s 


•33 


S| 




H 


3 


1 

1 

1 






CP 


1 

1 

CDI 








1 H| 


» 


» • 


• 


1 =n 


cj 


oi oo 


Oi 


| 








°i 


Cm 


<<s< 


i si 

1 i 








l 1 

1 01 








I H 








1 H 


cj 


U5 pq H uq 


1 h i 









co 

0) 

co a; co 

I I I I \ I M\ 
S CD £3 

CD 
Q 



• 

Oi 


9 

Oi 








9 9 9 


3 


3 


H 








Uh 




U 


Oi 

pq 


H 


CP 

2H 





P0P0 


rnrn 


j 




a a 


s a 


cnCP cp cp j 




W 


W 


HHHH 


1 


cncn 


cncn 1 


w w 1 




HH 


HH 


s | 

1 

1 

1 


9 


9 




1 

1 

1 




CJ 


COPM wm CJ > j 


CQ 


OhS 


SSpCfJOlpS 


04 


CP 


CP CP 


1 

OjCj Oj | 


H 


po 


(_5 u> * (H H3 H5 cj I 



! 



I 

I 



I 

I 



I 

I 

I 

I 



I 

I 

I 

I 

I 



I 



I 

I 



I 





0J 

Q) 

H 


1 

1 

I 










CD 




-Q 


1 


M 








CO 




fl* 


1 


m 








rrt 




•H 


I 


-p 








0 




M 


j 


3 










M 


3 


I 


3 3 




H 




M 


> 


> 


1 


=3 CP 




0 3 0 3 


3 


CD CP 






1 


M M3 


3 


OOOO H 


H 


O a -ph 


w 


a 


1 


OH 


3 


HHH-H 0) 


3 


H e CD 3 3 3 


H 


3 


0J| 


M CO 


M 


CO H CO -P 3 CD 


-p a 


4J 3 3 CD 3H 


cq 


M 


C0| 


0) M CO 




O 0 3 CP 


3 3 M 


CD a CD O > M -P 


3 


cn 


SI 


-3 HO 


H 


M3M3MCDC0 W 


000 


OH 3 M U CP | 3 


H 


0 


1 


-PH 3 M 


H 


*H 3-H 3 H O 


MH 


3 0 CD H pq CO oq 




M 


M| 


OO O 




03 H (T3 H H H ( D n H 


HH 


H M CO CD 3 CD 




CP 


°! 


3 M M M 


O 


O 3-3 CD O 


3 3 


M 3-3 OH 






i 


3H-H-P O M 


3 


qnqiri .h+j+j-mmq)3^-h 


3 3 CO 3 -PH 0 3 0 




3 


31 


CO 3 3H 0) 


*H 


-O^rOosaiM X U-H^H 0 0 


oq r — 1 -p cpO d'H 




•H 


3| 


H CD 3 -P 


a 


CH 0 :Mrt3 CO-H-H yj 'jxn 0 M 


3 3 3P3 MH M f 0 




3 


H| 


HHHH C0H 


M M a co aoc >i>h q uuid 


3 H O CO CD >*H O (DCN 3 CD 




2 


M) 


O UJ-H 3 


CD CD 


0 0 -P H PSd^P CO -p 


3 3 O^-CNH 3 CO 3 



fa | W Ul 3 -w-p U| O U| CJ CJ u O ro-r-i CJ UJ ol, M-P r-1 

CD) 3HH 0 3 CD O 3 -3-3 O M m -P-P-P3 3D3 

SI| H PP U H CDS? O-H 3 — I H 3 3 3-P-H 3 > S 

I B9^ W C7*3 P'OPHJ 3-3 3 3 -P CD S CP CD O OHH CD CD O Q 

CO I M3HC0 SOOU3<D3<UOJHO<D H M CD 0) O ^^flOOPU^CJTP 

-I cd m cu 3H 3 m-p o o co cd co co *h pcopoccp ocpcpo 33 

CD) -POC^CD-H H S CO S CO 3 3-3 -P-H-H > CQH <D 3 H -HSPCrfH 

HI O) O a CD M 3 3 3 3 O 3 3 3 Ml CO SJ 3 -P U-P-P-PH M 3 OH O 

H| 3M MH CO OH 3-H OPCJIPP CO HH CD a 3*H 03 3 3H CD H H 

tuj Mj(U (U U/T-l Url (OH H U( | PHO W tau M'H O' JO O' -H PHPHO 

•H| <J i— | u|-wh (J ro (0(0fU(00> (O H-H JHOI U ilj y U| PH ■* JHPHrHU t-H-H-H-n 

M| -M a) = M (tJ UJ=; > M rH rH M C3 M 3 (O MCO 0> 0) 3 -M -M -W -w U >fH r- 1 (J 

f0| -X O 0> O 0> W Ol/iawa ouotjic m M CDCQ-H >>,h-H-H cnrd nj >a03 

> t3+j 0) s»HrH MH'a s i-h«j t/iM a a) (D+j a tn rc to oi to d^m-x-x as a; 

| 0> rH-H r O'M<^VJ H U|i-t U| Ul-U O b C'M UH w c w c q q (0 y-w v y 

l c/vh occ Qi-c? oi tn o o o o ui o,g<m caom ukc ooM0>0)0iMcqc m-u 

I =0 OSSHQ-MB-HXMXMa Q«fl OH WU«3>MQrt;QQQQ[i|MWrt; 



I 



CD| 

HI 




25 H 




r - 


CN 


O 4 


fcH 


CD 




oq >-< 


r-CNCJ OOCP 


H| 


2C 


P? 0-1 


H 


2J 


atJ 


Oi 


UP 


UP 




vjt-i 


r^^UUlNUH 


3| 


M 


2 ^ Uh uc £_q 


O) 


ZD 


a) 


0-1 


H 


30-1 




c^iOO 


oq oq cd oq az pm h 


HI 




uq uq cq cd 


U4 


CD 


CD 


uq 


S 


z 0-1 


2; 


2 2 


2 2 3Pt-t tH tH 


HI 


CD 


HMM^ 


UP 


(JQ 


UP 


wx; 


O 


UP 


PO 


o-i PM 




3| 

>1 


«« 


*3 <5 


^=q 


^3 


•3 


3 CJ 


CJ 


CJCJ 




£-q Pq 


pq pq pq mm o-l o-i H 



I 

I 

I 

I 



I 

I 



I 

I 

I 



I 



I 



62 



r 






g 

o 













CO 
















CO 


CO 


H CO 


CO 


CO 


CO 


CO 


(01 






CD 


G 


rd rH 


H 


• — 1 


rH 


CD cn 


+J| 


CO 


CO 


a; 


<d 


-Q rd 


G 


rd 


rd CO 


DH 


•HI 1 1 1 1 


i i 1 i i \ i 


1 \ 


H 


•H 


O U 1 


U 


O 


CJ N 


H \ 


gi 


a 


a 




"a 


f — 1 c/j 


UJ 


Ui 


CO s 


OTO 


:g| 






cd 


fd 


*H H 




G 


G 


a) a 


1 






a 


« 


Oi 


Oi 


cm 




a 



CD 

H 

G 

H 

o 



cut 

OI 

H 

H| 



• • • « 

OJUJUJOJ 

* = d ^~C ^ 



SSS^(^H 






Uh 

OJ 



Pn 

LM 



U 



« 



Oi 

H 



O, 

ua^u v 



Oi 

H 



CM Op 
pU LJ 



ujuj :> 
<2: 3 



Oi 

WU 





U 


1 




CD 
















1 




-P 






CD 




G 






1 




cdH 






u 




G 1) 




CO 


1 


r— 


H O 






G 




rd Ha 


H 


0) 


1 


+ 3 


PMH 




CO 


CO 




a Hrd 


H 


rH 


1 


DD rd 


a-H 




u a 


CO 




0) G CD 


> 


-G 


1 


^G rd rd 3 0 


O '0 




rd 0 


CD 




CO < H 




(d 


1 


a a q-p-p g u 


O CD VJ 




H H 


u 




a g -w 


W 


•H 


1 


oaaannoH 


H H CD O CD 


CD 


O hj 


fH 




(dPQ g c/i 


H 


H 


a>! 


\G G d-HH H CO 


O rd -GHH 


H 


1 — 1 O-t 




e 


0) O PH 


C Q 


rd 


C0| 


CD CD CD OOHGd) 


-P d CO -PiG cn 


*T> 


•H O 


CD 


03 


Ul CD *H 0 


<J 


> 


:gi 


CJ OpC-WH 3 o 


O OMH-H rd G 


r-« 




-H 


CD 


■PG -P 


EH 




1 


O-i O fd 


rd UH G Pd G 


rd 


HJ 


-P 


H 


CO -P O 0) 




3 


HI 


O O O 0) CD COH 


V HH 




H CD H 




-P 


G cd a 




rd 


OI 


CO -P-P h 


H CD H G'H 


H 


•H 0 0 


H 


CO 


rsH ■ — 1 G 




SH 


1 


-P H H H cd rdt-CN G 


H CD H U H H O 


0 


03 -P 


O 


H 


0 OO'MH 




H 


Oi 


rd O O O CJ CJ 00 


O-U-H O rd 




-PMH H 




G 


H CD O 




o 


HI 


0)-H-H-H-H-H4J4-> >iH T3 a i-l+J 3 'U 


CO 


H H-H 


CO 




C0-PQ> 




H 


•HI 


z: Ul <S> V jrDTD Cfi01+J-P124!1H3^3 


CD 


•H O H 


'J 




CM CD G 




04 


HI 


CO CO CO a H-H-H H-H 


<a cy-^ w yj<-i 


a 


OOQ CH 


G 




G*H CD CJ 






ro| 


U CD CD CD-H-H OUUU 


O 't! -M -w CL) 3 fO 


H 


HH 


H 


-P 


-PH O-h’H 




a 


cd 


•HULIU OiUjUO-H a> CJ ecu > c/j -o 


CD O CD 03 


G f -D 


G id VP ' H 




•H 


s:| 


4H HUMLlpO O HH H'TJ H G 


CD > 


>iGP 3 


> 


rs-H O-H > -H 




rd 


1 


H X G CD 


U-iCJ rd HH-w 


CD 


GO — • 




- OH CJ -G p 




sc 


COI 


OD did) rd rd O > • H '4-1 O 0 3 


H m 


o3CJ CD rd GrvpOi dO ^ O 






» 1 


CD HH 


-p -P U <D a,T»G 


r A > 




G 


G ‘d WCV 






CD| 


U,q c C 1) Q) U U M !— 1 <— 1 OJ 0-3 G-i-’ <i)-H 


-P-P G 


•H-H-P 


pH CD -P CO 






HI 


co ca n3 ca ui ui <y m o-h 3 a oi \-m 0 


GGO-P r ao<d-PG r a g 






-HI 


ddHHG<2:H PO -HO 


H 


PH^^i H 3 


H -P O CD 






fd| 


HH'HM-IMH tq ta cq 


fdrH HHH 


H (J 


u»u;y< hu 


CD HCD CD ^3 U 






•HI 


UCIUUHWJD3M 


H ^ O O U fd <J‘ 


’H-H 


MM'C) PHH H H H-H CJ G 






HI 


CD CD G G G G 


03 > H-H 




G G HH 




O 3H 3 U3-P 0) 






rd | 


oa>cui)-t j 4J exj-p +j a> 3 


rd >ird W W 


>i(3 U) 


W >l3 P 






>1 


•H333C G-G-G-G-H O G-P G 3 0^0 3>td T>W V) U G 


3 dl U) G !fl 0) 



g g u cj cj *h upc/jhhh cdcducdcdvdcdodm cd ^ cd H'o cdh 
C do3rdo300rdo3rd <d-H c; Q; M aj (d a ^ a M M (/) H -P H U| H o3 CO <D 
«>>>cjos:s:s:s:^heh o ~ t > uaoco^^^ i 03 cm i cm < h m cp 



1 

CD| 
» — 1 1 
H| 


r-CNPO 

<3; 


COH 








00 

Up CJ H 




r— 

OJ 


CM CO 
OJ CJ 


H 

Hu 


rd| 




cp jj np op cj 




CD 


CD 


«d00 00 


OJ 


OJ 


CO H 


v > 


•HI 




2 ; cp cj cp ^ <*; 




H 


H 


WWW 


H 


U-J 


IV H 


V Up 


H| 


«4 *=U 


lj 


G) 


U4 


ra 


xw^ 


UP 


Up 


Up <q 


ujw 


fd| 

>| 


pCJCJC3 r O r J£3SZ;2UZi ^ 


2; 


CJ 


CD 


U 4 m U4 


H 


H 


HUp 


IWW 



63 



TABLE VIII 






i 


m 


in mm 


01 


















1 Wl 


CD 


c anq o g 


q q q 


C7i 
















1 +JI 


CD 


dH d CDH 


•H-H -13 


H 


01 


01 


C/1 


CQ C/J 


oi 


01 


V) 


| 'HI 


p 


•h xh P > 


> PH 


X 1 


X 


X 


X 


NX 


X 


X 1 


X 


i ci 


an 


neon oh 


i — i h n 


n 


a 


a 


a 


s a 


a 


a 


a 


i CD | 


CD 


d a d cd o 


O CD d 


















! I 

1 I 


Q 


CG kqu: 


h H C5 


















1 1 

1 d 
























1 ' — •! 




• I 


• 


• i 


• 


t 










• 


1 =>\ 




p i CQ 


o 


HP* 


o 


u 










o 


1 *01 


• 






















! °l 




cxzss 


SPS 


CQ3 S 


pq 


CD 


n 


S Z 


z 


z H 


pq 


1 si 

1 I 
























1 1 
1 <Dt 
























I aj 




CP Dui 


CP cu 


CP 


cp 


CP 


fp 


CP CP 


CP 


CP 


CP 


— >i 


« 


ppwpp 


H H H 


H H 


CP 


CP 


CP 


pp 


1—3 


H VJ 


H 



« 



« 

o 



HI 



H 

O 



>1 

H 

H 

O 



U 


! 


1 


n 




o o 








1 


1 


CD 




CD HH-HH H 




• 




1 


1 


tn CD 




G (1) CD -P 




rri 








d rp 




O >H>HH 




P 


0 


1 


i 


cd an 


a 


0 q G o 




o 


1 — 1 


1 


i 


P G 


d 


PH GH O 




o 


H 


1 


i 


c/i o d a 


CD 


H G O q OH 


O O H 


0 


d 


! 


i 


q o d 


p 


H d OPd HCD 


H H -H 


1 


*H 


1 


i 


d n cd o 


+J 


a a > 


O 


H P 


P 


CD] 


1 


•H G P 


01 


CD OH CDH 


r — 1 - — 1 O 


H O 


d 


c/> 




n 01HH 


q cd cd 


G 01 d 01 d o 


d d ' — I 


P OH 


> 


P>| 


dCNHH 0/ 


3 OH 


H 3 G H 


a a a o 


cncDH q 




| 


1 


h an 


O d H 


a ccicmqp 


P P d > 


qpp O 


a 


Pi 


1 0) 


h d d o 


nn 


pn o o 


o a o o 


H d o 


d 


o| 


1 H 


q q ancD 


pH 


0 ohhh-hh 


2(133 pH 


q 3H 


P 


1 


1 an 


HH CD pr- 


CN O O 


HH OH OH O 


CD HO 


•H O O 


cn 


oi q 


o gD 


on 


CD 0 O O 


Hp HOI 


a p h 


o 


G| <q 


CD HO 01 H 


H 01 


n d oi o c/i cu cq 


h -p h qn 


P O 01 o 


p 


•HI 


H G 


q CD CD 


CD O P O P CD 


HH 01 H 3 G 


O OP 


CP 


G 1 -H 


CTHHHH 


•H q 3 


o P GH GH =3 


-UO'^iOOai+J-HHG-H 




di u 


Q03UOU 


O OH UJ 


•m nnnnn 


•HOG O'O G G O-H-H O 


G 


CD! U CD H H HO 


no cd cd 


g d i d t cd 


OrH H Q O'O G G 


*H 


si 


1 H 


ZD CDH+-* ZD 


PH 


n o > x > hp n O OH O H HP 


p q 


d 


1 


1 oann =1 dl D4JHPP chh 


o cu 


h p q > p an 


mm q 




Oil 


tiU3 OH rd G G TiG5 O, 


V5^3 G>C cr>G OM-i O O O O 


o o and 




- 


1 cl) > OH >~H > 


C fdH 


G GH G-H q 


g>'-m g'-h h u^hh g a 



<U| 

HI 

m 

*1 

•HI 

Ml 

d| 

>1 

! 

I 

I 



n _c; on a cu 

H H 00 P» O pH 

g a a a 

MH WXHH U IN M 
U P CJHH U 03 q 
H H 01H PH 



0) OH cd 
p pn qj- 
a an ph 

-P PH q 
CO ro^H^ 
P P a 
O CD HdH 



o o*H o 
g q q o and nq na,d o d 

HdHcucDoadaappH 
d d CUP c /7 CD O O PH O 
> I HCOCOIOH I 



CQH H *HH OH OP (J H CD 

CD dn 01H 01H OS HOQ 01 n 01 

OH H -PH -PH Hd^H^H dH 3*HH G 
H Q P CJ P U cy q CO P UH UiH QW 
rQCUM.cU-MCUHnpCDCJOQ-M 03 «P 

>o q q UDGqqqqyj>a o 

d HO HO Hfd ^001003 OH P HH 
HH d Hd nd O 0 n n onH -P cd+J 

ojppqpqp o s O a cd p oi cd P 

OGPOPOPHP OH OH d H OH PH 
H^<q0<O<OHU-PU-PHfdH-P<qo 



I 



O! 

^1 


CO 


CN 


coi 


tH 


HI 




HH i- 




Mi 


d! 






Oi HH 


IH 


HI 


H 


WP»WWS 




fi 


PI 


JJ 


qMtrlPM 


HHH 


CJ 


d! 


c/J 


ujujua tn t-i 


HtriH 


t-l 


>1 











W 




r— 


CN 










OJ 2D 


SD 


o 


H 


H 




pm uq 


W 


m; 


(01 




H 


HO 


o 


<53 


PM 




PM 


PM Mi 


Mi 


tn 


H 


t-l D P> 


> 


P*P> 


P> 


p* 


pc 



64 



Main Program Variables (cont'd 
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Note: this legend applies tc all Tables listing progiaui variables. 
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Figure A. 8 Subroutine CHKINP Flowchart 
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Figure A. 9 Subroutine CHKINP Flowchart (cont'd. 
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figure A, 10 Subroutine CHKINP Flowchart (cont'd. 
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figure A„11 Subroutine DEFANG Flowchart 
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Figure A. 13 Subroutine WSHOCK Flowchart (cont'd 
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figure A. 14 Subroutine HATVEI flowchart 
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Figure A-15 Subroutine WATVEL Flcwcbart (coat'd. 
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AP|EN£IX_B 
SAflPLZ P BIN TO UTS 

This appendix presents copies ox the output from various 
runs cf the computer program listed in Appendix C. The 
varices saiiples illustrate the output from the following 
cpt i cis : 

(1) A Complete Print of a Water Pun 

(2) A Summary Print of a Water Bun 

(3) A Complete Print of an Air Run, and 

(4) A Summary Print cf an Air Bun 

(Note: in erder to reduce the volume of the tnesis / only a 

portion cf the output from the complete print cpticr.s is 
included . ) 
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SAMP IE 1 - COMPLETE PRINTOUT (WATER CAS 
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FNSIIT AT CCNE 5CRFACE = 1 243. 1082 KG/ M 3 

FAG CC EFFIC IE NT (CD) = 0.0630 
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SAMPLE 2 - SUMMARY PRINTOUT (RATER 
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CEEGA AT POINT 2 = 25.00C0 DEGREES 

STREAMLINE ANGLE AT POINT 2 = 4.1003 EEGR EES 
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E N I H AL £ Y AT CONE SURFACE = 1432619.0 J/KG 

DENSITY AT CONE SURFACE = 1328.2969 K3/M3 

DRAG COEFFICIENT (CD) = 0.1150 



SAMPLE 3 - COMPLETE PRINTOUT (AIR CASE) 
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ELSEBANN RADIUS A1 PCINI 3 = 745.5216 M/S 

X-CCCRD 1NATE OF ELSE MAN N CENTER = 1566.1 827 M/S 

X-CCCBEJNATE OF EUSEMAN N CENTER = 582.7557 M/S 
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SAKPIE 4 - SUMMARY PRINTOUT (AIR CAS I) 





















4 






























* 






























* 






























# 






























* 






























* 

’fr 




























































* 






























* 






























* 


• • 






























CO 






























H 




























* 


w 






















UJ 


W W 






w 






















N 


\\ 




* 


CO 






CO 


CO 0 














S 


ss 






M 






N 


MM 






















W 






s 
































CJ H 














r— 






* 


M 








C/3 M 














m 


or- 


CO 




O 






0 
















on 


CNO 


w 










CO^J 


H 














r- 


oo 


w 




X 






N- 








CO 








• 


• * 


w 


# 


W 






y-rsiin 


O 






W 








on 


Oin 


0 




<q 


CO 




• 0 


oon 






< 








ao 


oo 


w 


# 


S 


w 




XX 


in 






CJ 








co 


aocj 


w 




s 


M 


CO 




oa- 






CO 




s 


CO 


CO 






# 


0 


« 


N 


on in t 


0 • 










\ 


N 


\ 








CO 


0 


s 




ctin 






W 




CJ 


s 


=T S 




on 


* 




w 




r— 


rn^r 






25 


♦ 1 


w 




CN 




O 




< 5 ; 


Q 




t 


MX 






H cO 


w 


CO 


CO 


O 




0 








1— 




rn 






>0 CN 


W 


w 


*- N 


r^o 




0 




C/3 






r— 








W • 


w 


wo 


cj s 


• On 




• 


* 


H 


r— 


VO 


on 








WO fO 


3: 


WCN 


>o 


^in ii 


II II 


tops 






X 


O 










*qo 




ocn 


rn 






wo 


* 


O 


cn 


O # 






* 


• i 


vO 


CO 


wao 


• ct 


i X 


XX 


w 




25 


O 


OCO 


ll 






co 


<N 


2: 


a . 


aoo 


in Eh 


hh 


w 


* 


H 


1 


• O 






* 


s 


oo 


o 


r-* 


inn 


OH 


HH 


0 




3 


0 


cooo 


M 






o 


r-CN 


M 




cN^r 


on (j 


CJCJ 


w 11 


# 


O 


0 


M 


CJ 11 


II 


1 


W 


* CO 


H 


on 


* 


O 


oo 


w 




H 




W 








M 


< T ) m 


M 


CO 


r- 


• — i 


HH 


0 




M 




H 


P-j M 


M 


Q3 


u 


i> w 


W 


(N 


CN 


M 


MW 






0 




O 


II HCJ 


II CJ 


H 


W 


rn qz 


25 


o 


in 


>■ 


>P> 


OH 


* 


M 




W 


03 








CJ n 


CJ 


• 


T— 


II 




« 25 






11 


W II 


M CO M 


M M 




CO 


w 


O 


r- 




w 




OH 


* 


1 






U 03 


CJ w 


1 




II Q « 




CN 


ll 


CN O 


S2; 


2TO 






w 


w 


<wo 








w 


W 




II 


II 


^=q< 


H 


* 


Q 


H 


own 


W^CO 


wco 




w 


W CO 


*=c 


II 


X 


H H 


ss 






W 


O 


tOQ 


wo 


03 




w 


« ll os 


H 




XH 


25 CN 25 


MM 


H 


*«• 


W 


S 


in s x 


O CJ M 


OM 


W 


M 


W #0 


m 


II 


CO I — I 


M M 




ll 




CJ 




^ron 


u j 25 


U j M 


M 


i-3* 


M HU ) W 


M 




H W 


W M 25 


M M 










^M 


MW 


W 






<< w^r 


21 


WH 


CJO 


H 250 


MW 


O W 




M 


X 


CJ 


W-=CCJ 


M CJ 


co 


CO 


WO W 


H 


W25 


OH 


H H 




H 


* 


« 


W 


32 O 


2 : 


2: 


H 


w 


wco CJ 




OH 


HW 


HOS 


(M W 


HO 






H 


CJ H 


OQH 


OH 




23 


WW II <q 


H 


25 W 


M5> 


«4MO 


CJCJ 


2: 25 


■» 


25 


U5 


II <W 


U2< 


w *** 


25 


W 


25 w s 


W 




> 


CJ 




H<53 




W 


W 


S> 


O 




o 




WWW 


H 




w 


WH 


EhH 


O 




W 


> 


W 


HOW 


HW 


w 


> 


H WW£3 


< 


2:H 


SH 


w «< W 


2 : 2 : 


HW 




PQ 


1 


HS^j 


rt-CO W 










W 




*=c<q 


w < 


WM 


2 : 






H 


CJ ^ ^ 


W 


3D 


V j 


H 


20S2:m 


o 


M 


M 


S X M 


252 : 


MM 




C/3 


S 


25 M M 


XMS 


MM 


M 


O 




u 


MX 


wo 


o MM 




<c w 




< 


M 


<«W 


M WO 


H<3 


32 


W 


W CM M 


h 


CJt-1 


H O 


•MM 25 


HH 






33 


CO 


MH 


H 25 


OH 


M 


25 


^UnXUj 




M M 


MS 


CJ M 


Ml 


<q 


•* 






WUjUj 


CJ M 


U J M 




M 




CJ 


HUj 


MM 


SOW 


VJCJ 


OM 




M 


M 


CJ M W 


WM O 


COH 


*■ 




UJ l/J WM 




M 25 


M ^ 


UHM 


CJCJ 


*M W 




25 


25 


CJ M M 


HM W 


MS 










M M 


W^ 


«=q m<s 


1 1 


SM 




CJ 


CJ 


M H W 


MM 


W M 


■* 




UDujm 




MW 


IMS 


S P> M 


X X 


c JUJ 




CJ 


CJ 


w 3M w 


j 2: 


WM 



8A 



EFNSJ3Y AT CCNE SURFACE = 6.2334 K3/M3 

EFAG CCFFFICIENT (CD) = 0.7549 
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PBCGBAh IIS1ING 

Ihis appendix certains a complete listing on the fullp 
documented rain program named CGNEF10W and all or its func- 
tions ard suhreutires. Ihe functions and subroutines 
included in this appendix are as follows; 

( 1) function DTCB 

(1) function BTOI 

(E) function DZEMG 

(*J) function SOSAIE 

(5) Subroutine CEFINP 

(6) Subroutine HSECCK 

(7) Subroutine KA3VEI 

(8) Subroutine EP^ 
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CALCULATE FINAL AIR CONDITIONS AT THE CONE SURFACE 
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EE5FCNSE IS CCEBECT. TAKE ACTION DIfiEClED 
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NCTE; THE FUNCTION HAS A FEATURE WHICH TESTS TC EEI ERMINE IF THE DE- 
FLECTIOR ANGLE HAS BECOME LARGER THAN 1 RAEIAN. IF IT HAS, THE 
MAIN PROGRAM'S EXECUTION IS STOPPED DUE TC A FATAL ERROR- IN 
ACCORD A KCE WITH KINNEY AND GRAHAM'S, EXPLOSIVE SHOCKS IN AIR, 
THE DEFLECTION ANGLE CAN NEVER EE LARGER THAN 1 RADIAN- 
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E£GIN (FUNCIICN SCSAIR) 

SOSAIR = LSQRT(G1 * (VMAX**2 - SV£L**2)) 
RETURN 
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